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SUMMARY 


One -dimensional, steady-state, compressible, viscous flow relations 
are presented which permit the determination of flow conditions at any 
radial position in a ducted helicopter blade. The relations are required 
for estimating the performance of proposed helicopter jet-propulsion 
systems which involve ducting air or gases through the blade from root 
to tip. A limited number of calculations over a wide range of helicopter 
operating conditions and relative duct sizes are also presented. The 
"choking" problem in the straight duct is discussed. 


INTRODUCTION 


The use of jet engines as propulsive units for helicopters is being ' 
investigated by many agencies. Several different power-plant configura- 
tions are under consideration, including those which require the ducting 
of air or gases through the length of the helicopter blade. Analytical 
relations are presented herein for evaluating the conditions of the ducted 
flow at any point of its travel from the blade root to the tip. Such 
information is required in order to determine limiting quantities of flow, 
optimum duct sizes, overall losses, and flow properties (i.e., static 
pressure, stagnation pressure, and Mach number) of the gases delivered 
to an engine unit at the blade tip. 

In subsonic ducted flow the action of wall friction and turbulence 
losses in producing decreasing density and increasing Mach number is well 
known and has been described adequately in the literature (e.g., ref. l). 
The limit to this action is the attainment of a Mach number of -unity or 
the "choking" condition. The choking problem may require careful con- 
sideration in the design of ducted helicopter blades because the passage 
length relative to the passage diameter necessarily will be large and the 
entering Mach number may be in the high subsonic range. 

The effect on the ducted flow of the centrifugal compression due to 
the blade rotation will be to raise the density and lower the Mach number. 
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and thus counteract the effect of friction. The net effect of the two 
opposed actions will depend on the particular combination of flow con- 
ditions, duct geometry, , and helicopter operating conditions and must be 
evaluated for each combination. 

The analysis presented herein is limited to steady-state, one- 
dimensional, compressible, viscous flow and applies only to straight 
ducting. 
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SYMBOLS 

p 

duct cross-sectional area, ft^ 
speed of sound, fps 

2 Tg 
T t 

specific heat at constant pressure, Btu/lb/°P 

specific heat at constant volume, Btu/lb/°F 

lift 

duct hydraulic diameter, 

Perimeter 

friction factor from reference 2 
acceleration due to gravity, ft/sec 2 
Joule’s constant, TJd ft-lb/Btu 
duct-flow Mach number 
rotor tip Mach number, AR/ a o 

duct-flow stagnation pressure, lb/ft 2 abs 
duct-flow static pressure, lb/ft 2 abs 

radial distance from center of hub to blade element, ft 
blade radius, ft 

Reynolds number of ducted flow at station 1, Pi^i^l/Pl 


rotation parameter. 
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t a (M) 

tb(M) 

^c(M) 

t d (M) 


static temperature, °R 

duct-flow stagnation temperature, °R 

duct-flow velocity, fps 

ratio of blade-element radius to rotor -blade radius, r/R 
ratio of specific beats, Cp/Cv 

mass density, lb-sec 2 /ft i4 ' 
rotor angular velocity, radians/sec 
coefficient of viscosity, lb-sec /ft 2 
Mach number function for heat term 

Mach number function for friction term 

Mach number function for centrifugal-force term 

Mach number function for area-change term 


Subscripts: 

0 atmospheric 

1 duct station at hub center 

2 duct station at rotor tip 

x any radial station 


ANALYSIS 


For steady-state, one -dimensional, compressible, viscous, ducted 
flow subjected to a centrifugal force in the flow direction, as in flow 
through the length of a helicopter blade, the rate of change of momentum 
of a mass segment of the flow may be set equal to the algebraic sum of 
the forces on the segment (Newton's second law of motion) as follows: 
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pVA dV = -A dp - kk | V dr + i^pAr dr - p dA 


( 1 ) 


Equation (l) assumes a constant molecular weight and specific heat, no 
drag-producing bodies in the stream, and no changes in mass flow rate. 

On the right-hand side of the equation, the first term represents the 
pressure force, the second, the friction force, the third, the centrif- 
ugal force due to the duct or blade rotation, and the fourth, the pres- 
sure force due to a change in area. 

Equation (l) may be converted into terms of stagnation temperature 
and Mach number through use of the following: the continuity relation] 

the perfect-gas law; and the expression relating stagnation temperature, 
static temperature , and Mach number. 

The converted equation is as follows : 


^ M (l + 5«2)(l + lriM2) aTt + mil M 2 ) 

2T t (l - M 2 ) a^ + D (l - M 2 ) 


M^l + 7 ~ 1 M 2 ) 2 ^ n(l + M 2 ) ^ 

6 (cp - c T )jT t (l - M 2 ) 1 - M 2 ^ 

A study of reference 3 will show that equation (2) is equivalent to equa- 
tion [19] of that reference except for terms which the previously enu- 
merated assumptions eliminate and the centrifugal-force term which was 
not considered in reference 3- 

Equation (2) may be nondimensionalized and expressed in terms of 
a rotation parameter B as follows: 
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where 


ta(M) 






1 - M 2 


(5b) 


(3c) 


* c (m) = 


7MI1 + 



l - M 2 



(3d) 


b = % a (|) 2 (^) (3e) 


and 

*b(M) 

27 M 2 

The solution of equation (3a) requires the use of additional rela- 
tions such as expressions for the duct-flow stagnation temperature and 
the duct cross-sectional area in terms of radial position in the duct r 
or the differential equations thereof. Such an expression for the stagna- 
tion temperature would permit the evaluation of the first term on the 
right-hand side of equation (3a), which may he classified as a heat term. 
The exact nature of such an expression depends on many design details and 
on flow and operating conditions, and its detailed analysis is beyond the 
scope of this paper. A few general remarks, however, are in order. 

The influence of the heat term would be to increase the Mach number 
if heat were added and to decrease the Mach number if heat were transferred 
from the ducted flow. Heat would be added to the flow if the blade duct 
were used as a combustion chamber, in which case the relation of stagna- 
tion temperature to duct length would be a function principally of the 
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space rate of combustion and associated parameters . A certain amount of 
heat transfer from the flow will generally occur with or without combus- 
tion because the ducted-flow temperature will be higher than the local 
stagnation-temperature recovery on the external blade surface if it is 
assumed that the ducted flow is compressed mechanically before it enters 
the blade duct. At any station in the duct the amount of heat transfer, 
and thus the temperature of the ducted flow, would depend on the heat- 
transfer coefficient, the stagnation-temperature recovery of the ducted 
flow, and the stagnation-temperature recovery on the external surface of 
the blade. The heat-transfer coefficient would be determined to a large 
extent by the blade and blade-duct design. 

The second term on the right-hand side of equation (3a) accounts for 
wall friction effects and would tend to increase the Mach number. The 
friction factor f is a function of Reynolds number and the degree of 

roughness of the duct surface. In the Reynolds number range of 10^ to 

107, the data of reference 2 indicate that the following empirical rela- 
tions apply: 

For commercial pipes (steel, cast iron). 


f » 


0.0247 


%e 


0.141 


For smooth surface (glass, copper, drawn tubing). 


f = 0.0236 

w °' 155 
^Re 


These relations indicate that at a given Reynolds number the friction 
factor for rough surfaces is about 30 to 40 percent greater than for 

smooth surfaces. Increasing the Reynolds number from 10^ to 10? produces 
about the same percentage decrease in the friction factor. It is shown 
in a subsequent section that changes of these magnitudes in friction- 
factor values do not affect the solution of equation (3a) appreciably. 

The third term of equation (3a) accounts for the influence of the 
blade rotation, which is to compress the gas and reduce the Mach number. 
The term contains a rotation parameter B, defined by equation (3e), which 
involves the rotor tip Mach number, the ratio of duct hydraulic diameter 
to blade radius, and the ratio of atmospheric temperature to the duct- 
flow stagnation temperature. It will be noted that the definition of 
rotor tip Mach number, 
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depends only on the tip speed due to rotation and the atmospheric speed 
of sound. 

The fourth term in equation (3a) accounts for the effects on Mach 
number of changes in cross-sectional area and is dete rmin ed by the duct 
geometry. 

Equations (3a) to (3f) provide a means for obtaining the Mach num- 
ber variation along the blade duct. The relations used in converting 
equation (l) to equation (2) may be combined to provide an expression 
for the static-pressure ratio across a length of duct: 
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(5) 


Equation ( 5 ) may be converted to an equation for stagnation-pressure 
ratio through use of the relation between static pressure, stagnation 
pressure, and Mach number to produce 
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(6) 


It is of interest to note that for adiabatic flow in a constant- 
area duct not subjected to centrifugal forces, the first, third, and 
fourth terms on the right-hand side of equation (3a) are eliminated. 
The resulting equation can be integrated to give 


7 - 1 2 
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( 7 ) 
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Equation ( 7 ) represents an end condition for the helicopter hlade-duct 
performance in which the rotational speed has been reduced to zero, and 
is of use in obtaining an overall view of the effect of the blade rotation. 


NUMERICAL CALCULATIONS 


In order to obtain a better understanding and a graphic view of the 
solution expressed in equation ( 3 a.) , a few numerical evaluations were 
made which were directed toward investigating the choking problem. The 
calculations were limited for simplicity to the case of adiabatic flow 
in a constant-area duct. The limitation to adiabatic flow is justified 
on the basis that heat transfer from the flow would tend to decrease the 
Mach number and thereby prevent choking, so that the adiabatic-flow case 
is conservative in this respect, and heat addition in the blade duct 
(combustion) is not believed to be a strong possibility. 

The calculations were based on assumed values of the rotation param- 
eter B (see eq. (3e)) of 0, 0.0002, 0.0007, and 0.0014. The value of 0.0014 
corresponds to a ratio of blade radius to duct diameter R/D of 26.7 for 
a rotor tip Mach number Mr of unity and a temperature ratio Tq/T^ of 

unity. These conditions are representative of an estimated maximum prac- 
tical value of B or an estimated min i mnm practical value of r/d. A 
zero value of B corresponds to the other extreme of no rotational speed 
or infinite blade radius relative to the duct diameter. For adiabatic 
flow and no rotational speed, equation ( 7 ) is applicable. 

For each assumed value of B, equation (3a) was evaluated for adia- 
batic flow for duct-inlet Mach numbers (at the hub center) of 0. 3, 0.6, 

and 0.8 at an inlet Reynolds number of 10'. The value of the friction 
factor f was taken from reference 2 and corresponded to commerci al 
cast-iron pipe. The rough-surface friction factor was chosen in order 
to obtain a conservative answer relative to choking. 

The adiabatic, constant-duct-area version of equation (3s) was 
evaluated numerically by a step-by-step integration process to obtain 
finite changes in Mach number for a series of values of A(r/D). The 
values of the Mach number functions of equations (3c) and (3d) were 
taken from precalculated curves presented in figure 1. A curve of equa- 
tion ( 3 b) ha 6 also been included in figure 1. In the Mach number range 
from 0.4 to 1.0 the reciprocal of the Mach number functions is presented 
to avoid having the ordinate approach infinity near a Mach n umb er of 1.0. 

The changes in Mach number in the blade duct were determined in 
increments of length equal to 4 duct diameters until a total duct length 
equal to 32 diameters was reached. The increment of length equal to 
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4 diameters was found to be small enough so that the Mach number could 
be calculated -with an error less than 0.5 percent- In order to obtain 
these accuracies, however, it was necessary to base the calculation of 
the Mach number change across a given increment on the Mach number at 
the center of the increment as predicted from the change in Mach number 
across the preceding increment. After completi n g the calculation for a 
given increment it was necessary, for the first 8 diameters of duct 
length, to repeat the calculation on the basis of a corrected estimate 
of Mach number at the center of the increment. 


RESULTS 


The variation of duct-flow Mach number, as obtained by these proce- 
dures, is presented in figure 2. For the cases corresponding to the higher 
values of rotation parameter (B = 0.0007 and 0.0014) the centrifugal effects 
more than compensated for the friction effects, so that substantial decreases 
in Mach number occurred through the length of the duct. This result was 
obtained even at an inlet Mach number of 0.8, for which the stationary- 
blade curve (B = 0) indicates choking or a Mach number of 1.0 at a duct 
length of 7*8 diameters. For an inlet Mach number of 0.8, a value of B 
of 0.0002 produced a Mach number of 1.0 at a length of 10. 5 diameters; 
however, a value of B of 0.0003 did not produce a choking condition 
and resulted in a net decrease in Mach number over the blade length. 

The curves of figure 2 specify the Mach number variation in the 
blade duct for given values of the parameter B, inlet Mach number, and 
inlet Reynolds number. These in l et conditions, however, would actually 
be determined by comrpressor operating point, nozzle-exit area, and other 
factors. Therefore, in order to change the inlet conditions or move 
from one group of curves to another in figure 2, a change in one or more 
of these factors is implied. 

The influence of Reynolds number or friction factor on the Mach num- 
ber was investigated briefly. Calculations were made for B = 0.0007 at 

a Reynolds number of 10 6 in order to compare with results for the value 

of 10? of figure 2. Reducing the Reynolds number to 10^ is equivalent 
to increasing the friction factor about 40 percent and resulted in a 
3-percent increase in the Mach number. This relatively small influence 
of Reynolds number justifies ignoring small changes in Reynolds number 
in calculating the flow along the blade duct. 

The Mach number at a duct length of 30 diameters as a function of 
inlet Mach number is presented in figure 3, which is a cross plot of 
figure 2. The effect of the rotation parameter B on the duct exit Mach 
number is clearly illustrated. An approximate choke line intersecting 
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each B curve has been drawn in figure 5 hy fairing through points corre- 
sponding to choking at some location in the duct. The line was faired 
through a point for B = 0 calculated from equation (j) , through a point 
for B = 0.0005 determined by cross plots of figure 2, and asymptotic 
to the vertical at an inlet Mach number of 1.0, according to the reasoning 
that as the rotation parameter approaches large values the inlet Mach num- 
ber must approach 1.0 in order to obtain choking. The choke line repre- 
sents a locus of the minimum values of inlet Mach number for which a Mach 
number of 1.0 can be obtained at some location within the 50 -diameter 
duct length. As the rotation parameter is increased from a value of 
zero, the radial location at which a Mach number of 1.0 will occur moves 
from the 50 -diameter location upstream toward the hub, and the minimum 
inlet Mach number increases from a value of O .657 and approaches 1.0. 

If the duct length were increased, the minimi mi inlet Mach number for 
choking when B = 0 would decrease according to equation ( 7 )- 

The stagnation-pressure and static-pressure ratios over the 50-diameter 
length are presented in figure 4, which was obtained through use of equa- 
tions ( 5 ) and (6) and the Mach number values of figure 5- The intercepts 
of the pres sure -ratio curves on the ordinate axis are indicative of the 
pressure difference over the blade length due to centrifugal force with 
no internal air flow and were calculated from the expression 


yB 

^2 = ^2 = e 2(D/R) 2 
Ptq Pi 


which was obtained by integrating analytically between stations 1 and 2 
the expression for the centrifugal force on an increment of gas mass. 
With values of rotation parameter B of 0.0007 and 0 . 0014 , stagnation- 
pressure ratios of 1*55 and 2.41 were obtained with no flow, and the 
stagnation-pressure ratio decreased slightly with increasing inlet Mach 
number because of increasing friction losses. The increasing divergence 
of the curves for static-pressure ratio and stagnation-pressure ratio 
with Increasing inlet Mach number is a result of the increasing differ- 
ence between inlet and exit Mach numbers. (See fig. 5*) 


CONCLUDING REMARKS 


The differential equation presented for steady-state, one-dimensional, 
compressible, viscous flow was used in a step-by-step integration process 
to determine the Mach number changes for adiabatic flow in a ducted heli- 
copter blade over a wide range of helicopter operating conditions and 
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relative duct sizes. With adiatatic flow in a constant-area duct 30 diam- 
eters in length, a Mach number of 1.0 will not he attained for inlet Mach 
numbers of less than O. 657 , which corresponds to the case of no blade 
rotation. As the duct length increases, this value of minimi mi inlet Mach 
number for choking will decrease. For a given duct length, as the rota- 
tion parameter increases, the minimum inlet Mach number for choking 
increases, and the choking point moves toward the inlet end of the duct. 

The maximum pressure ratio across a length of duct occurs with no 
flow. The stagnation -pres sure ratio decreases with increasing inlet 
Mach number because of increasing friction losses. For the range of 
rotation parameters assumed, a maximum stagnation-pressure ratio of 2.41 
across a 30 -diameter duct length was calculated for the adiabatic -flow 
case. 


A change by a factor of 10 in the duct-flow Reynolds number, which 
determines the friction factor and friction losses, produced only a 
3-percent change in the calculated duct-flow Mach number for the adiabatic- 
flow case with a moderate value of the blade rotation parameter. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 5, 1953- 
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(a) M = 0 to 0.6. 


Figure 1.- Mach number functions defined by equations (3b), (3c), and (3d). 
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(b) M = OA to 1.0. 


Figure 1.- Concluded. 
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Duct length , r/D 

Figure 2.- Relation of duct-flow Mach number to duct length for constant 
values of B = l^ 2 (D/R) 2 (To/Tt) • % e = 10?. 
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Figure 3 .- Net Mach number change for duct length of 30 diameters. 
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Duct- f/ou Mach number at r/D-O 


Figure k.- Stagnation- and static-pressure ratios for duct length 

of 30 diameters. 
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